Measurement of Magnetic-Field Structures in a Laser- Wakefield Accelerator 
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Experimental measurements of magnetic fields generated in the cavity of a self-injecting laser- 
wakefield accelerator are presented. Faraday rotation is used to determine the existence of multi- 
megagauss fields, constrained to a transverse dimension comparable to the plasma wavelength ^ Xp 
and several Xp longitudinally. The fields are generated rapidly and move with the driving laser. 
In our experiment, the appearance of the magnetic fields is correlated to the production of rela- 
tivistic electrons, indicating that they are inherently tied to the growth and wavebreaking of the 
nonlinear plasma wave. This evolution is confirmed by numerical simulations, showing that these 
measurements provide insight into the wakefield evolution with high spatial and temporal resolution. 
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The last few years have witnessed tremendous progress 
in the field of laser-driven electron acceleration. Electron 
beams with monoenergetic spectra can now be produced 
with high-power laser systems [1|, potentially providing 
a source of brilliant secondary radiation y. Further 
improvements have now been made in terms of energy 
[3| and stability [J]. In these experiments, the high- 
intensity laser pulse drives a plasma wave ^], which 
evolves into a 'bubble'-like structure with dimensions of 
order Xp — 2ttc/llIp both transversely and longitudinally 
[if. The 'bubble' compresses the laser pulse both spa- 
tially [6( and temporally [7[ provided that it exceeds the 
power for relativistic self-focusing P > Per = 17{ncr/ne) 
GW. This can lead to an increase in laser intensity even- 
tually leading to wavebreaking and injection of electrons 
into the wave's electric field. 

This high charge of accelerated electrons (> 100 pC), of 
ultrashort bunch duration {rb < t; ~ 50 fs), thus consti- 
tutes an extremely high current (> kA). Since the current 
is confined to a dimension much smaller than the plasma 
wave (r < /^rn), it will generate a large associated az- 
imuthal magnetic field, B^ = jiqI /2T:r. The field is not 
neutralized over distances comparable to the magnetic 
skin depth (c/up) and over associated timescales (l/wp), 
which are larger than those of the electron bunch. Hence 
these fields will serve to collimate the beam. 

This is not the only magnetic field associated with the 
wakefield. In one dimension the coherent motion of elec- 
trons in the plasma wave do not produce a magnetic field, 
as the current due to the plasma oscillation jc — nev ex- 
actly cancels the displacement current produced by the 
resulting charge separation jd — eo ^ • However in three 
dimension, electrons are pushed ofF-axis and return via 
the sheath of the 'bubble', so that there is a circulation 



of current in the wakefield which generates an azimuthal 
magnetic field [8J. In the 'bubble' regime the field can be 
calculated by considering the displacement current as- 
sociated with the moving electron void B^ ~ ^at^^^ 

^ f 2 a{x — ct) 

[9|, where rb is the blowout radius of the void and Ex 
is the longitudinal electric field [10|. Hence the size of 
the field grows rapidly with plasma wave amplitude, and 
thus with increasing laser intensity. Thus measurements 
of these fields could provide important information which 
is difficult to obtain from highly non-linear wakefield ac- 
celerators [5] , except from their radiation properties |ll| . 

In this Letter, we present the first experimental evi- 
dence for the generation of multi-megagauss azimuthal 
magnetic fields in a laser wakefield. The evolution of 
the magnetic field indicates the rapid growth and sub- 
sequent breaking and decay of a large amplitude non- 
linear plasma wave. Simulations confirm this evolution 
and show that the field is a combination of that due to 
the displacement current of the plasma wave as well as 
that associated with the current of relativistic electrons. 

The experiments were performed with the 10-TW 
Ti:Sapphire JETI laser at the lOQ, Jena. Pulses of 
wavelength A — 800 nm, energy E^ = 800 mJ, and 
duration r = 85 fs were focused by an //6 off-axis 
parabola to intensity / ~ 3 x 10^* Wcm~^ onto the ris- 
ing edge of an interferometrically characterized helium 
gas jet from a 1-mm diameter nozzle (Fig.[T]). Accel- 
erated electrons were either detected by a fluorescent 
screen shielded by a 30 /xm Al-foil to block laser light 
and a CCD-camera or alternatively by a high-resolution 
magnetic spectrometer. Quasi-monoenergetic features 
peaked around 50MeV were observed in the electron 
spectrum as in [12!|. A synchronized probe pulse was gen- 
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FIG. 1: Experimental set-up. An electron beam is gener- 
ated by a self-injecting laser wakefield accelerator. A linearly 
polarized probe pulse traverses the interaction region and ex- 
periences Faraday rotation. The interaction region is imaged 
onto two CCD-cameras, each analysed with a polarizer. 



erated from the transmission of the laser pulse through a 
1:100 beam splitter. The probe was telescoped down and 
guided into the interaction chamber via a variable delay 
13[. Due to group- velocity dispersion when prop- 



line 

agating through glass, the probe duration was slightly 
increased to 100 fs. The interaction region was backlit 
by the probe, which was then split by a non-polarizing 
beam splitter and imaged by a high-quality //lO lens 
with a magnification of 10 onto a pair of CCD-cameras, 
equipped with interference filters. The spatial resolu- 
tion was 10 /zm in the transverse and 30 /im in the laser 
direction, determined by the imaging system and due 
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FIG. 2: Two polarograms in false colour for same laser shot 
with polarizer angles: a) Spou = -1-5.9°; b) Q^ovi = —4.1°); 
c) distribution of 0rot deduced from the intensity ratio of a) 
and b). Most pronounced rotation occurs in a small region of 
(55±5) /imx(35±5) ^m indicated by black arrows; d) line-out 
between the red arrows shown in c); e) associated (jirot. 



to the motion blur caused by the orthogonal motion of 
the main beam relative to the probe. In addition, po- 
larizers (extinction ratios of (1 — /3i) = 4 x 10^"^ and 
(1 — /32) = 4 X lO"**) were positioned in front of each 



141 was used to 



camera. A Nomarskii interferometer 
measure the electron density. 

The i3,^-field distribution during the interaction was 
measured using Faraday rotation [15|. When B^p has 
components parallel to the propagation of the probe, the 
probe polarization is rotated by <i>mt, (Fig. [I]). The ac- 
cumulated rotation is (/)rot ~ -^ / ''^c{f)Bip{r) ■ ds 

taken along the path i of each ray of the probe through 
the plasma. Here, rid- = eQiUcUiij. / e^ the critical den- 
sity for the probe light of frequency Wpj. and wavenumber 
/cpi- = Wpr/c. To measure the distribution of (/)rot, the 
two cameras were carefully aligned to take two images 
of the same interaction region. The two polarizers were 
detuned in opposite directions by 0poi,i and 0poi,2 from 
the positions crossed with respect to the initial probe po- 
larization. The intensity transmitted through polarizer i 
(where i = 1, 2) can then be described by Mains' law. 



/pou = /o [1 - ft sin^ [90° + e, 



pol,i 



t]] 



(1) 



where Iq is the initial probe intensity. Taking two simul- 
taneous polarograms eliminated shot-to-shot fluctuations 
of the probe-beam profile including plasma refraction ef- 
fects, thus increasing the sensitivity of our measurements. 

Fig. [2^) and b) show two such polarograms. The 
main pulse has entered the plasma from the left; inten- 
sity differences can be seen in the centre of the images. 
This is even more pronounced in the distribution of 0rot 
(Fig. [21;), deduced from the ratio of the intensities of the 
two polarograms using eq. ([1]). The regions where the 
strongest rotation occurs are confined to the center of 
the image, being symmetric around the laser axis and 
having a longitudinal extent (55±5) /Ltm and a transverse 
diameter of (35±5) /zm. The main pulse was situated 
slightly in front of this region, indicated by the dotted 
ellipse in Fig.[5J;). The position of the main pulse was 
determined by the appearance of ring-like interference in 
the images on separate shots where the probe intensity 
was reduced to a level comparable to the side-scattered 
light from the main pulse. As the interference can only 
be generated when the probe and scattered light overlap 
in space and time, the center of these rings marks the 
instantaneous position of the main pulse. A lineout of 
the intensity ratio through this region and of the associ- 
ated Faraday-rotation angle are shown in Fig.[2jl) and e) . 
The appearance of the intensity variations was correlated 
with the generation of an electron beam. 

Two-dimensional (2D) simulations of the interaction 
were performed using the particle-in-cell code OSIRIS [l6( . 
The plasma density profile was set to be gaussian (1/e^ 
radius of 550 /im) with a peak density of 4 x lO^^cm"'^ = 
0.023ncr as in the experiment. A laser pulse of gaussian 



temporal profile with r = 70 fs (FWHM), and vacuum 
focus of 6 /im (1/e^ intensity radius) corresponding to 
/ = 5.4 X lO^^Wcm^^ was set to focus 430 ^m before 
the peak density. A slightly higher initial intensity was 
chosen since 2D simulations tend to underestimate effects 
such as self-focusing which are important in the experi- 
ment. The simulations were performed in a box size of 
800 X SOOc/wo with 4000x1000 grid cells with 2 macro- 
particles per cell. The box moved at the speed of light in 
the lab frame to allow higher resolution. 

Fig. [3^) shows the electron density and the laser inten- 
sity from the simulation at the time when B^ is max- 
imum. Due to the evolution of the laser, the plasma 
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FIG. 3: Simulation results: a) electron density (blue) over- 
laid with laser intensity (red/yellow), b) longitudinal current 
density, c) B,p d) 92 /xmx92 /xm section of Fig.[2t), e) intensity 
modulation obtained from the numerical simulation including 
all effects affecting the image formation in the experiment, f) 
vertical line-outs along the dotted lines in d) and e). 



wake has evolved from an initially multi-period wave 
to a 'double-bubble' structure. Self-modulation of the 
pulse has occurred in both the longitudinal and trans- 
verse dimensions; the resulting pulselets contain suffi- 
cient energy to drive a highly non-linear wake structure 
12l |l7[. Both 'bubbles' have a strong B,^-fLe\d and are 
of a similar size, Fig. [3]:). The first 'bubble' has a field of 



magnitude B^p ~15MG, whereas for the second 'bubble', 
B^ ^35 MG. However, only the second 'bubble' has an 
on-axis current of trapped electrons generating a i3;^-field 
(Fig.[5)D) indicating that there must be additional contri- 
butions to the field which is present in both 'bubbles'. 
This is evidently due to the displacement current. Both 
B^ contributions are orientated in the same way (left 
handed in direction of laser propagation). Even though 
the peak B^ lies within the 'bubble', where the density is 
low, in general there is not complete cavitation and the 
Btp spreads into the high-density walls of the 'bubble' 
leading to a measureable Faraday-rotation signal. Note, 
however, that the B^-fields are mainly confined to the 
'bubble' structure due to return currents carried by back- 
ground electrons flowing around both 'bubbles'. 

For a quantitative comparison between experiment and 
simulation the numerical results shown in Figs.|3^)-c) 
were post-processed with a routine explicitly taking into 
account experimental effects influencing the image for- 
mation. Based on a ray-tracing routine accounting for 
refraction of the probe in the plasma, we further included 
the limited spatial resolution of the imaging system to- 
gether with the limited temporal resolution due to the 
flnite probe pulse duration of 100 fs and the blurring 
caused by the perpendicular motion of the main pulse 
within the duration of the probe. These effects serve to 
increase the apparent size of the feature and also reduc- 
ing the signal's amplitude. A comparison between the 
S^-field feature's spatial extent and a vertical lineout of 
the intensity distribution from the experiment and the 
post-processed simulation data are shown in Figs.[3Jl)-f). 
A good agreement is found both in size and amplitude. 
For a field of multi-megagauss strength to be generated 
by the accelerated electrons, implies a current of such 
electrons > 10 kA. The simulation shows that such large 
value of field and current are only possible because the 
pulse duration of the electrons is much shorter than the 
plasma wave duration (of order ~fs). 

To separate the contribution to B^ by the wake dis- 
placement current from that from the electron-current, 
the individual quasi-particle momenta in the simulation 
was recorded. A particle was considered to be trapped if 
its forward momentum px > 7 MeV/c, corresponding to 
the velocity of the plasma 'bubble' for rio > 9 x 10^* cm~^ 
and Aq = 800 nm [9J. From the current density of all 
trapped particles the associated azimuthal magnetic field 
is calculated. It increases rapidly at the same time as 
the total B^ confirming that for our experimental con- 
ditions the generation of B^p is strongly correlated with 
the trapped electron current. However, the accelerated 
electron alone only accounts for approximately half of 
the total i?;p-field (as shown by the circles in Fig.|4]), con- 
firming the importance of the B^-iie\d due to the 'bubble' 
itself. 

By controlling the delay between main and probe 
pulse, the temporal evolution of the magnetic fields 



was visualized (Fig.|4^-f). Starting from the first im- 
age where no clear signature can be seen, the feature 
appears and reaches a maximum rapidly, but then be- 
comes weaker again for later times. The distance be- 
tween the 'bubble'-positions in two subsequent images 
which were delayed by 200 fs stays approximately con- 
stant at (60 ± 5) fXTa leading to a propagation velocity of 
the 'bubble' of (3.00 ± 0.25) x lO^ms^i, matching the 
laser's group velocity in the plasma. 

In Fig. 13^), the evolution of the peak Faraday- rotation 
angle, which is proportional to the B^ and which has 
been measured at different times in the plasma, is com- 
pared to the peak B^ from the simulation. After a rapid 
increase close to the peak of the density, both quantities 
show a slower decrease towards the end of the plasma. 
This rapid onset and the subsequent evolution both of 
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FIG. 4: a)-f) Time evolution of B^. 
270 fim) shows distribution of (^rot- 



Each image (560 /im x 
g) Peak (/)rot (squares) 



as function of position compared to corresponding B^ from 
simulation (red). Relative plasma density profile is indicated 
by dashed line. Circles denote B(p-contribution from trapped 
electrons. Error determined by difference between peak value 
of (irot above and below laser axis. 



the rotation and the B^-fie\d imply a rapid steepening 
and breaking of the plasma wave and subsequent injec- 
tion of background electrons into the wake. From these 
measurements, it is evident that the wave breaks only 
once, rather than a continual injection which is seen in 
simulations at higher intensity. The experimental mea- 
surements decay faster than B^p in the simulation. Hence 
it is apparent that in the experiment either the confine- 
ment of accelerated electrons and/or the plasma wave 



amplitude decays much faster than in the idealized con- 
ditions of the simulation. This is the first visualization of 
such non-ideal evolution of a plasma wave, and can help 
to explain the often large discrepancies between experi- 
ments and simulations. 

In conclusion, we have measured the evolution of mag- 
netic field structures during the interaction of a high- 
intensity laser pulse with underdense plasma. This is 
closely related to the process of non-linear plasma wave 
formation and its breaking into a 'bubble'-like structure. 
The sudden onset of the i3<^-field is a clear signature of 
the non-linear steepening of the plasma wave and, for 
many experimental conditions, a signature of the subse- 
quent injection of electrons into the wakefield. By reduc- 
ing the plasma density below the threshold for electron 
self-injection our technique could also be used to detect 
the B— fields of the bubble alone. Our diagnostic gives 
experimentalists a powerful tool that can easily be imple- 
mented into existing set-ups to provide information of the 
acceleration process with high resolution. The results can 
directly be compared with numerical simulations which 
presently are the basis of our understanding of the ac- 
celeration processes in the plasma, eventually enabling 
researchers to identify crucial experimental parameters 
required for a stable interaction. This will help in op- 
timizing laser-driven electron acceleration, bringing the 
prospect of table-top sources of brilliant and ultra-short 
XUV and x-rays closer to reality. 
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